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The angle-resolved photoemission spectrum of the organic conductor TTF-TCNQ exhibits an 
unusual transfer of spectral weight over a wide energy range for temperatures 60-?^ < T < 260K. In 
order to investigate the origin of this finding, here we report numerical results on the single-particle 
spectral weight A{k,u!) for the one- dimensional (ID) Hubbard model and, in addition, for the ID 
extended Hubbard and the ID Hubbard-Holstein models. Comparisons with the photoemission data 
suggest that the ID Hubbard model is not sufBcient for explaining the unusual T dependence, and 
the long-range part of the Coulomb repulsion also needs to be included. 

PACS numbers: 71.10.Fd, 79.60.Fr, 71.20.Rv, 72.15.Nj 



The low-dimensional interacting systems receive at- 
tention because of their unusual electronic proper- 
ties Q. In this respect, the high-resolution Angle- 
Resolved Photoemission Spectroscopy (ARPES) mea- 
surements on the quasi-one-dimensional organic conduc- 
tor tetrathiofulvalene-tetracyanoquinodimethan (TTF- 
TCNQ) have provided evidence for non-Fermi liquid be- 
havior in this compound H, Q IS- In particular, the 
ARPES experiments have found that the single-particle 
spectral weight at the Fermi wavevector A:f is transferred 
over an energy range of ~ 1.3eV of the Fermi level in 
the TCNQ-derived band, as the temperature T decreases 
from 260K to 60K HQ. In a Fermi liquid the spectral- 
weight transfer would have occurred within ~ fc^T of 
the Fermi level. Here, we investigate the origin of this 
unusual ARPES data and its meaning for the electronic 
structure of TTF-TCNQ by using the Dynamical Den- 
sity Matrix Renormalization Group (DDMRG), Quan- 
tum Monte Carlo (QMC) and the exact diagonalization 
methods. 

There are various possibilities as to what might be 
the origin of the anomalous T dependence of the single- 
particle spectral weight at the Fermi level, A{kp,uj), 
in TTF-TCNQ: (i) It has been suggested that the T- 
dependence of the ARPES data can be explained within 
the one-dimensional (ID) Hubbard model [^Q. In this 
case, the anomalous T-dependence of A{kF,uj) over the 
conduction bandwidth has been attributed to the strong- 
correlation effects. Indeed, by using the Bethe-Ansatz 
solution, the photoemission spectrum has been fitted ex- 
cellently to the dispersion of the spinon and holon bands 
of the ID Hubbard model with the parameters t = OAeV 
for the hopping matrix element and U — 2eV for the 
Coulomb repulsion. The recent observation of the 3kp 
structures in A{k, to) by the ARPES also supports this 
picture, (ii) An alternative point of view is that an ex- 
tended Hubbard model with long-range Coulomb repul- 
sion is necessary, particularly because the screening of the 



long-range Coulomb repulsion is expected to be weaker 
for the surface layer of TTF-TCNQ. (iii) Another possi- 
bility is that the electron-phonon interaction, in addition 
to the strong Coulomb repulsion, plays a role in produc- 
ing the unusual T dependence. In this paper, our goal is 
to differentiate among these possibilities. For this pur- 
pose, we present DDMRG and finite-temperature QMC 
results on A{kp,uj) of the ID Hubbard model. In addi- 
tion, we present exact-diagonalization results for the ID 
extended Hubbard model which includes a near-neighbor 
repulsion V and DDMRG results for the ID Hubbard- 
Holstein model. 

In the following, we show that, above a characteris- 
tic temperature determined by the effective magnetic ex- 
change Jeff, spectral- weight transfer takes place over a 
wide energy range in the ID Hubbard model. This is 
similar to the T dependence observed in the ARPES ex- 
periments. However, below this temperature, the weight 
transfer is negligible. We find that the Hubbard param- 
eters t = OAeV and U = 2eV give too large a value for 
Jef f , and with these parameters it is not possible to ex- 
plain the ARPES T dependence. On the other hand, the 
T — exact-diagonalization results on the ID extended 
Hubbard model show that the nearest-neighbor Coulomb 
repulsion increases the bandwidth of the spinon and the 
holon excitations. This can lead to a smaller value for t 
for fitting the ARPES dispersions, and a reduced value 
for Jeff, hence, giving better agreement with the ARPES 
data. In addition, the DDMRG results on the Hubbard- 
Holstein model show that, at T = 0, the electron-phonon 
interaction influences A{kF,uj) only at |w| ^ 1.3eV for 
physical values of the phonon frequency loq. However, it 
remains to be seen how the electron-phonon interaction 
influences A{k,uj) at finite T. The main finding of this 
paper is that the ID Hubbard model is not sufficient for 
explaining simultaneously the T dependence and the dis- 
persion of the photoemission spectrum of TTF-TCNQ. 
We suggest that it is necessary to include at least the 
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long-range part of the Coulomb repulsion. 
The Hubbard hamiltonian Hq is defined by 



Hn = -t 



i.<7 



{4,aCi+i,a + h.c.) + (1) 



where Ci^a {c\ c) annihilates (creates) an electron with 

spin a at lattice site i, Ui — n^^j + n^^i, rii^cr — cl^Ci^a, 
t is the hopping integral, and U is the on-site Coulomb 
repulsion. We will consider the case of electron-filling 
{n) = 0.60, since the fiUing of the TCNQ band is 0.59. We 
note that A{k, w) of the ID Hubbard model was studied 
with the QMC '0,'7','i,|3 and the DDMRG as weU 

as with the Bethe ansatz 

Within the DDMRG method, the single-particle spec- 
tral weight is obtained at T = from 



A(k,uj) = Im( 



H 
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-Cfc,T 



(2) 



where Cfe_| annihilates an electron with wavevector k and 
spin tj |0) and i?o are the ground state and the eigenen- 
ergy, respectively, and 7 is a small positive number. The 
DDMRG results were obtained with the open boundary 
conditions. At finite temperatures, we obtain A(k,uj) 
by using the determinantal QMC technique This 
method yields the single-particle Green's function along 
the Matsubara time axis, from which A{k, lo) is obtained 
by the maximum-entropy analytic continuation 0|. We 
have checked the convergence of the maximum-entropy 
results as the statistics of the QMC data improved. In 
the following, Aik^uj) will be plotted in units of t^^ . In 
order to determine the characteristic temperature of the 
magnetic correlations, we present QMC data on the uni- 
form magnetic susceptibility x{l ~^ 0) where 



Jo 
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(3) 



with m| = n^.i — n.i^x. 

We also present exact-diagonalization results on 
A{k,uj) for the ID extended Hubbard model. 



Ho 



(4) 



where V is the nearest-neighbor Coulomb repulsion. The 
A{k, oj) of the ID extended Hubbard model was previ- 
ously studied with the exact diagonalization technique 
[1^. In addition, we present DDMRG results on A{kp,uj) 
of the Hubbard-Holstein model defined by 



Hhh ^Hq + ujoY^ ^1^* + 9^(4 + bi)ni 



(5) 



where bj (bi) is the creation (annihilation) operator for 
a dispersionless phonon at site i, ujq is the phonon fre- 
quency, and g is the electron-phonon coupling constant. 
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FIG. 1; (color online) Temperature dependence of the single- 
particle spectral weight at the Fermi wavevector, A{kF,uj), 
for the ID Hubbard model with (a) [/ = 5t and (b) U = 8t. 
In (c), the change in the ARPES intensity with respect to 
T = mK, AI{kF,uj), is plotted (in arbitrary units) for TTF- 
TCNQ (reproduced from Refs. y, 3 with permission) . 



We first discuss DDMRG and QMC results on the ID 
Hubbard model. Figures 1(a) and (b) show the T de- 
pendence of A{kF,uj) for U = 5t and 8t. The DDMRG 
results were obtained for a 60-site chain with 36 elec- 
trons, while the QMC results are for a 32-site ring with 
(n) — 0.60. In the DDMRG calculations we have used 
a finite energy broadening 7 = 0.05t. For U = 5t, we 
observe that the spectral-weight transfer occurs over an 
energy range w 2t of the Fermi level, as T decreases from 
0.5i down to 0.125t. However, the weight transfer is neg- 
ligible between T = 0.125t and T = 0. For t w AOOmeV 
glil, T = 0.125i corresponds to « 600K. Hence, in 
this case, the amount of spectral-weight transfer between 
T = GOOi^T and T = OK is negligible, which disagrees with 
the ARPES results. In order to study the dependence on 
U /t, in Fig. 1(b) we show results for U — 8t. Comparison 
of Figs. 1(a) and (b) shows that the transfer of weight at 
low T is enhanced for U = 8t with respect to U — 5t. In 
Fig. 1(c) we show the change in the ARPES intensity at 
kp, AI{kF,uj), for TTF-TCNQ HQ. These results rep- 
resent the photoemission intensity arising mainly from 



3 




0.1 0.2 0.3 0.4 0.5 -1 -0.5 

T/t u (eV) 



FIG. 2: (color online) Temperature dependence of the uniform 
magnetic susceptibility xi'i ~^ 0), plotted in units of t^^ , for 
the ID Hubbard model. We use these data to determine the 
effective magnetic exchange. 



the TCNQ-derived band. Here, the solid curve repre- 
sents the difference in I{kF,(^) between 260K and 60K. 
The transfer of intensity is also observable as T is low- 
ered from lAOK to 6QK (dashed curve). It is considered 
that, at low energies l^^^l^ OAeV, the interchain hop- 
ping becomes important Hence, here we will dis- 
cuss the energy range — l.SeT^ ^ w ^ —OAeV. In Figs. 
l(a)-(c), we observe important differences between the 
ARPES data and the numerical results. With U — 5t, 
it is not possible to explain the low temperature scale of 
the weight transfer. In addition, we observe that it is 
not possible to explain the energy range. In the ARPES 
data, the intensity at —1.3eV ^ w < —OAeV decreases 
with decreasing T. In the ID Hubbard model and for 
t = OAeV, the decrease of the spectral weight occurs 
for —O.SeV < w ^ 0, while the holon peak develops at 
u! w — 0.85eT^. Hence, it is not possible to explain the 
T and u! dependence of the ARPES data using the ID 
Hubbard model. 

In order to determine the effective magnetic exchange 
for the ID Hubbard model, in Fig. 2 we show the 
T dependence of the uniform magnetic susceptibility 
x{q 0). We obtain a value for Je// by making compar- 
isons with the ID Heisenberg model. In the ID Heisen- 
berg model, as T is lowered, xiQ ~^ 0) starts to decrease 
with the development of the 2kF magnetic correlations, 
and the maximum of xil 0) occurs at w 0.64Je// 
(isf. Obtaining Tm from Fig. 2, we find Jeff ~ 0.35t and 
« 0.26t for U = 5t and 8t, respectively. Hence, we have 
Jeff « leOOif for t ^ OAeV and U = 2eV. The results 
displayed in Fig. 1(a) for U = 5t show that the transfer of 
weight over the wide energy range occurs for T ^ Jeff/3, 
implying that this process depends on the development 
of the short-range magnetic correlations. Apparently, the 
long-range magnetic correlations are not playing a role 
in this case. At this point, we also note that xil ~^ 0) 
of bulk TTF-TCNQ is strongly suppressed by the fluc- 
tuations of the Charge-Density- Wave (CDW) gap as T 



FIG. 3: (color online) Exact-diagonalization results on the 
single-particle spectral weight A(k,uj) at fc = and k — n/i 
for the ID extended Hubbard model which has onsite and 
nearest-neighbor Coulomb repulsions U and V, respectively. 
These calculations were performed for U — 2eV on a 16-site 
ring with 10 electrons corresponding to (n) — 0.625. Here, 
spectra obtained for t — 0.25ey and V = 0.5eV are compared 
with that for t = 0.4ey and V = 0. 



decreases from the room temperature down to the CDW 
transition temperature Tcdw = SSif [lillr^ . 

The density-functional-theory calculations deduce that 
the hopping parameter t — 0A75eV for bulk TTF- 
TCNQ, while the analysis of the ARPES data yields the 
Hubbard parameters t = OAeV and U — 5t This en- 
hancement of t for the surface layer has been attributed 
to a possible tilting of the TCNQ and TTF molecules 
at the surface. In this paper, we have seen that the pa- 
rameters t = OAeV and U = 2eV give Jeff ~ 1600i^, 
which is too high to explain the T dependence oi A{k,uj). 
At this point, we suggest that the long-range part of 
the Coulomb repulsion might play an important role. 
For demonstration, we present exact-diagonalization re- 
sults on A{k,uj) for the ID extended Hubbard model 
with U — 2eV. Figure 3 compares A{k,uj) obtained for 
t = 0.25eV and V = 0.5eV with that for t = OAeV 
and V" = at wavevectors fc = and fc — tt/A. For 
the 16-site lattice, fc = 7r/4 is the closest wavevector to 
kp. We observe that, at fc = and for t = OAeV and 
V ^ 0; the holon and spinon branches are located at 
~ — 0.68ey and « —0.27eV, respectively. For parameters 
t — 0.25eV and V = O.beV, these structures are located 
at similar energies. Hence, it is possible to reproduce 
the locations of the spinon and holon branches by using 
a reduced value for t within the ID extended Hubbard 
model. This behavior is also observed at fc = 7r/4. Figure 
3 also shows that V induces incoherent spectral weight 
at higher |cj|. For the ID extended Hubbard model, we 
expect Jeff oc 4i^/(C/ — V). These results suggest that 
taking into account the long-range part of the Coulomb 
repulsion in fitting the ARPES dispersion can lead to a 
reduced t and, hence, might reduce Jeff and the charac- 
teristic temperature for the single-particle weight trans- 
fer. However, it is necessary to calculate A{k,uj) for the 
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FIG. 4: (color online) DDMRG results on A{kF,uj) of the 
Hubbard-Holstein model for uio = 0.2t and electron-phonon 
couplings 3 = and g = 0.2t. These results were obtained 
for U = 5t and broadening 7 = 0. It on a 20-site chain with 
12 electrons, for which kp ~ 6n/21. This figure shows that, 
at T = 0, the electron-phonon interaction with ujo = 0.2t 
infiuences A{kF,(^) for < 1.2t. 

ID extended Hubbard model at finite T. 

In order to study the effects of the phonons within 
the presence of the Coulomb interaction, we next present 
DDMRG results on A{kp,uj) for the Hubbard-Holstein 
model. In TTF-TCNQ, inelastic neutron scattering ex- 
periments 18] have revealed longitudinal acoustic and 
optical phonon modes. The optical branch is weakly 
dispersive and has a frequency of about lOmeV. Here, 
we present results for U = 5t and an Einstein phonon 
mode with ujq — Q.2t and g — 0.2t. In Fig. 4, we 
show A{kp « 67r/21, oj) for these parameters on a 20-site 
chain with 12 electrons. For comparison, we also show 
results for 5 = 0. For the 20-site chain, the finite-size 
gap Aps near fc^? is « 0.3t, and in this figure we have 
shifted the spectrum by Aps/2 so that the main peak in 
A{kp,uj) occurs at cj = 0. Because of the finite broaden- 
ing 7 — OAt, we do not observe the opening of the CDW 
gap. However, we see that the spectral weight at w w 
is transferred to ~1.2t ^uj ^ ~ 0.2t. We have performed 
similar calculations using luq — O.bt, where multi-phonon 
peaks are induced at energies « —ujq, — 2wo: ~3wo; etc. 
with respect to the location of the main peak. For phys- 
ical values of ujq « lOmeV, we expect the changes in 
A{kp,uj) at T = to occur at <C 1.3eV. These results 
suggest that the electron-phonon coupling is not partic- 
ularly important for investigating the spectral weight lo- 
cated at —1.3eV < a; < —OAeV. However, it is still 
necessary to study A{k,uj) of the ID Hubbard-Holstein 
model at finite T, since the electron-phonon interaction 
can influence the T dependence of the magnetic correla- 
tions and, hence, the T dependence of A{k,uj). 

In summary, we have studied of the ID Hub- 

bard model in order to investigate the unusual T depen- 
dence of the photoemission intensity of TTF-TCNQ. We 
have also presented T = results for the ID extended 
Hubbard and the ID Hubbard-Holstein models. We find 



that in the ID Hubbard model the transfer of the single- 
particle spectral weight takes place over a wide energy 
range above a characteristic temperature, which is too 
high to explain the ARPES data. We have shown that 
the long-range part of the Coulomb repulsion can lead to 
a reduced value for J^jf and, hence, to a better agree- 
ment with the ARPES data. Our results on the Hubbard- 
Holstein model show that, at T = 0, the electron-phonon 
interaction does not influence A{kp, lu) over the wide en- 
ergy range observed by the ARPES. In conclusion, these 
calculations give theoretical support to the notion that 
the anomalous T dependence of the photoemission spec- 
trum of TTF-TCNQ is due to the strong-correlation ef- 
fects as suggested by Claessen et al. However, we 
also emphasize that the ID Hubbard model is not suf- 
ficient for explaining the unusual T dependence, and at 
least the long-range part of the Coulomb repulsion needs 
to be included. 
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